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In this talk we construct BPS-saturated domain walls in supersymmetric QCD, for any values of the masses 
of the chiral matter superfields. We compare our results to those already obtained in the literature and we also 
discuss their range of applicability, as well as future directions that would be desirable to explore in order to 
achieve a complete understanding of supersymmetric gluodynamics as a step in improving our knowledge of how 
QCD works. 
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1. Introduction 

Supersymmetric (SUSY) gluodynamics, the 
theory of gluons and gluinos with gauge group 
SU(A^c), shares a lot of common features with 
ordinary QCD, with the advantage of being su- 
persymmetric. Therefore it seems a good idea to 
exploit the supersymmetric properties of SUSY 
gluodynamics in order to try to understand bet- 
ter how QCD works. Although this is an idea 
which has been around for many years, the re- 
cent and positive results of Dvali and Shifman [Q , 
which we shall briefly explain in detail, triggered 
a great deal of activity in these past three years. 

Our starting point will be the Lagrangian for 
SUSY gluodynamics, which is given by 
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where go is the gauge coupling constant, the 
usual gluon field strength tensor and A is the 
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gaugino field. This theory has an axial U(l) sym- 
metry which is broken down to Z2Nc, the chiral 
symmetry, by the anomaly. Moreover, the theory 
is strongly coupled: at a scale A gaugino con- 
densates (which, from now on, we shall denote as 
(AA)) will form. This will break the chiral sym- 
metry from Z2Nc to Z2 and, clearly, the gaugino 
condensate is the order parameter associated with 
the symmetry breakdown. 

The value of the gaugino condensate was cal- 
culated by various methods a few years ago [^. 
The Nc degenerate vacua of our theory are given 

by 



(AA) = (TrA"A„) = A'^e 
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The presence of degenerate vacua leads immedi- 
ately to the formation of domain walls interpo- 
lating between them. The key result of Dvali and 
Shifman [|l| was to realize that if these walls are 
of a certain kind, i.e. BPS-saturated, then they 
preserve one half of the supersymmetry and their 
energy density is exactly calculable and given by 

e = ^ |(TrA''Aa)oo - (TrA'^A„)_oc^| , (3) 

i.e. the difference between the values of the con- 
densate at the two vacua -situated at spatial 
infinities- between which the wall is interpolating. 
This is very interesting because it has been sug- 
gested p that BPS-saturated domain walls would 
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play an important role in the D-brane description 
of N=l supersymmetric QCD (SQCD). Also, re- 
cent work claims that, in the large Nc limit, these 
walls are BPS-saturated states but so far it 
has not been fully proved for finite Nc configura- 
tions. In any case the question to be answered, 
before making any further connections to other 
interesting topics, is what the nature of these do- 
main walls is. 

In order to do that we are going to discuss 
first of all different effective approaches to SUSY 
gluodynamics and the problems associated with 
them. Then we shall study in detail the Tay- 
lor, Veneziano, Yankielowicz (TVY) approach |^ 
which is the one used in our calculation. Finally 
we present results and draw a few conclusions. 

2. Towards an efTective theory of SUSY 
gluodynamics 

A first and very complete description of SUSY 
gluodynamics was provided by Veneziano and 
Yankielowicz Q (VY), in terms of a composite 
chiral superfield S whose lowest component is the 
gaugino condensate 
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The Lagrangian for this model is given by the 
usual supersymmetric structure, 
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with /C the Kahler potential and W the superpo- 
tential. Throughout this talk we shall be work- 
ing with canonically normalized fields, therefore 
JC = (5S')i/3. The superpotential for the VY 
model is given by 
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The structure of W is uniquely fixed by the 
anomaly and the symmetries of the theory. Again 
A is the scale parameter. It has been widely 
shown that the VY superpotential describes very 
well the vacuum structure of SUSY gluodynam- 
ics, however this is not the case when we try to 
describe the dynamics of the theory and, in par- 
ticular, when we try to construct domain walls. 



It was pointed out not so long ago [[7| that 
this Lagrangian is not explicitly Z2Nc invariant. 
Moreover, the scalar potential derived in global 
supersymmetry, i.e. V{S) = JC~^\dW/dS\'^ is 
multivalued due to the presence of the logarithm 
in W. That is, a physical state and its equivalent, 
just rotated by 27r, would give rise to completely 
different values of V, which is totally unaccept- 
able. 

In order to cure this problem, several solutions 
have been proposed. In Refs |7|j|] the idea of glued 
potential was introduced and developed. Essen- 
tially it consists of adding a Lagrange multiplier 
to the Lagrangian in such a way that the scalar 
potential is divided into Nc sectors which, glued 
together, result in a single-valued Z2Ar^-invariant 
theory. The problem with such an amendment 
is that cusps will inevitably form at the join- 
ing point of each sector with its nearest neigh- 
bours. Any configuration interpolating between 
two vacua will necessarily cross a cusp and it is 
doubtful whether it would be possible to correctly 
interpret the energy density associated with it. 
Therefore we conclude that the VY model, even 
when properly modified, is unable to provide us 
with a good description of domain walls in SUSY 
gluodynamics. 

Other solutions have been proposed to deal 
with the problems associated to the presence of 
the logarithm in the VY superpotential. For ex- 
ample, in Ref. 0| it was proposed that the origin 
of such problems was related to leaving behind 
some relevant degrees of freedom when deriving 
the VY efective Lagrangian. It is then suggested 
that introducing new degrees of freedom in the 
model, in particular a glueball order parameter, 
would be enough to enable us to construct well 
behaved domain walls. This is certainly a pos- 
sibility worth considering, however we will not 
be following this approach but that of Taylor, 
Veneziano and Yankielowicz ^ (TVY) , which we 
describe in detail in the next section. 

3. Our model: the TVY approach 

Given all the problems of the VY model when 
trying to construct domain wall solutions of the 
equations of motion, it is convenient to work with 
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a slightly more complicated model which respects 
all the symmetries of the theory. This consists of 
adding Nf pairs of chiral superfields, Qt to the 
VY model. Then, below the condensation scale 
we shall have matter condensates, Mj = Q^Qj, 
as well as gaugino condensates. The Lagrangian, 
given again by Eq. (||), has the following super- 
potential 



From now on we shall work in a flavour-diagonal 
basis, i.e. m-^ = Slnij. Our goal will be to con- 
struct domain walls between the vacua of this 
model, and then try to obtain the limit of large 
masses where the theory should tend to SUSY 
gluodynamics. In other words, we start off by 
working in the weak (Higgs) phase of the model, 
and we will try to extrapolate the results obtained 
to strong coupling. 

Therefore the first step consists of defining the 
vacua of the model. These are given by 
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where, and from now on, we are taking A = 1. 
That is, the vacuum values for the gaugino and 
matter condensates are aligned (* denotes values 
at the vacuum). Also, one can easily evaluate the 
superpotential to find out that it is proportional 
to the vacuum values of the fields, i.e. VF* = 
-(2/3)7V,^,. 

A very important point which we are going to 
discuss in detail is the paths the fields take when 
going from one vacuum to another. As mentioned 
before, we are using the matter fields to restore 
the lost Z2Na invariance of the TV Lagrangian 
and ensure that, when building domain walls, we 
do not cross the logarithmic branch. In order to 
formulate this in a more precise way, let us sup- 
pose that we are building a wall between vacuum 
a and vacuum 6. Then we define the field trajec- 
tories as follows 



M:\h - e^(*-2^"-)7WlU ,z = l,. 
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where LUi, the windings of the matter fields, have 
to be integers, in order to fulfill the alignment 
condition at the vacua. To ensure that the log- 
arithmic branch is not crossed one must impose 
the following condition, coming from the fact that 
the phases in the logarithm of Eq. (Q) cancel out 



iN,-Nf)S + yiS-2iTUj,)^0 
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As an immediate result it follows that S = 
2Trk/Nc where we define k = ^i- Note that, 

if all the windings are equal to one, then k = Nf. 
This is an important fact to which we will return. 

Let us now briefly discuss the method used to 
calculate the domain wall profiles. We assume 
that the walls spread along the xy plane, therefore 
the profiles are calculated along z. This is done 
by minimizing the energy functional which, for 
a generic set of fields, , with Kahler metric 
ai = K^x-x, > looks hke 

eafc = / dzigl^,X'^aJ + {9ir'^^^') , (11) 

where Wi — dW/dX^. This equation can be 
rewritten as 

eab = Re(e*^(Wb- Wa)) 
+ 2 



dz{gi)-'{gid,X''~e'''W^) 



(12) 



where e'^ is an arbitrary phase. In fact, if we 
choose it in a clever enough way, i.e. e"*'*'"'' = 
{Wb -Wa)/\yVb-Wa\, then it is easy to see that 
we can put a lower bound on the energy density 
of the wall that interpolates between vacua a and 
b, which is 



eab > \Wb~Wa\ 



(13) 
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This is the so-called BPS bound; if the bound is 
saturated the wall to which it corresponds is a 
BPS-saturated domain wall. An immediate con- 
sequence of the BPS condition being fulfilled is 
that the first term in Eq. ( |l^ ) becomes zero, in 
other words, the equations of motion become first 
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order, which represents a significant simplifica- 
tion from the numerical point of view. Further 
details of this can be found in Ref . . 

In particular, for the TVY model we are work- 
ing with, the BPS equations are 

,5W 



K-ssdzS 



95* 



dW 



(14) 



In order to simplify the analysis, from now on 
we will work with a Kahler potential which is 
canonical for the dimension one matter fields. 
For the case of degenerate windings, which we 
will analyze next, this is equivalent to taking 
Tr(MM)i/2 modulo factors of ^/Wj. The case of 
non-degenerate windings is slightly more involved 
and will be fully discussed in the future . Also 
we choose ^ = —\{5 + tt) = — ^ where, we 
stress again, k is the sum of all the windings as- 
sociated with the matter fields. From now on, 
in order to find solutions to these equations, we 
adopt the following parametrization 



S{z) = \SAR{z)e'^^'^ 
Ml{z) - |A4|p,(z)e^"'(-) 



(15) 



where R, (3, pi and at are real functions of the 
coordinate z. Note that the moduli of the fields, 
R and pi, are normalized to one at the vacua. The 
boundary conditions, when going from vacuum j 
to vacuum j + 1, are 



5 — > 5*6 



i2TT-, 



;27r(- 



(16) 



We also have a constraint, which is a consequence 
of the equations of motion, Eqs (pH), namely 



Im(e*'^>V(5,M^)) = constant 



(17) 



which, evaluated at the centre of the domain wall 
(i.e. z — 0) becomes 



Rf) 



1 - 



N 



Nr. 



(lni?o-l) + ^^lnp, 

i 

^-Y^[-ir^p,,^cos(n^) . (18) 



This constraint will turn out to be extremely rel- 
evant to the understanding of the results we are 
presenting in the next section. 

4. Results 

4.1. Equal windings 

To start the discussion of the results, we are 
going to present those corresponding to models 
where the windings uji are all equal to one (re- 
member that, in such case, k = Nf). This is the 
standard choice that other authors in the liter- 
ature have made, and it is therefore the one we 
should consider in order to compare with previ- 
ous results. In this case the constraint given by 
Eq. (na) becomes 
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Let us summarize the results. The case Nf = 
Nc — I was analyzed by Smilga and Veselov in a 
series of very interesting papers |l0|] , and recently 
reanalyzed by Binosi and ter Veldhuis [0. They 
found that the TVY model had BPS-saturated 
domain walls up to a certain value of the mass 
of the chiral matter condensate, m^. Between 
TO* and another, higher, value of the mass, to**, 
they found domain wall solutions which were not 
BPS-saturated. Finally, above to** there were no 
solutions at all. The values of m*, m** seem to 
depend on the value of Nc and they decrease when 
the latter increases. Moreover there seem to be 
two branches of BPS-saturated solutions which 
merge at to*. These results are very well illus- 
trated by Fig. 1 of the third paper in Ref. pC| ]. 

A couple of years ago, two of us analyzed 
the case of Nf — 1 flavour. We found that there 
were BPS-saturated solutions for any value of the 
mass parameter m. This can be seen in Fig. 1 
where we plot the magnitude of the gaugino con- 
densate, R, versus the spatial coordinate, z, for 
several values of m. As the mass increases it can 
be seen that the profiles tend to a unique one, rep- 
resented by the thick line. We shall return to that 
point later on. These results were subsequently 
confirmed by the work of Ref. flit] . In general it 



Figure 1. 

Plot of the magnitude of the gaugino condensate, R, as a 
function of the spatial coordinate, 2, for Nc = 3, Nf = 1 
and different values of the mass of the matter condensate: 
m = 2 (dotted), m = 20 (dash-dotted), m = 100 (dashed), 
m = 200 (solid). 



Figure 2. 

Plot of the phase of the logarithm of the TVY superpo- 
tential (see Eq. (^), in units of 2n, as a function of the 
spatial coordinate, 2, for the TVY model with Nc = 5 and 
Nf = l,m = 1000 (dashed) and Nf = 2,m = 200 (solid). 



is possible to make the following statement: 

• for Nf/Nc > 1/2 there are BPS-saturated 
domain walls only up to to* 

• for Nf/Nc < 1/2 there are BPS-saturated 
domain walls for any value of m 

Also, a more thorough study of the Nf/Nc < 1/2 
case [|l^ leads us to claim that, for Nf > 1, these 
domain walls will cross the logarithmic branch 
and therefore we are uncertain of their physical 
meaning. We illustrate this point in Fig. 2, where 
we plot the phase of the logarithm that appears 
in the second equation of Eqs (0) as a function 
of the spatial coordinate, z, for a TVY model 
with Nc = 5. As we can see, for Nf — 1 and 
a large enough mass, the total variation of the 
phase is never bigger than 27r. However, for Nf — 
2 and m = 200 the overall variation of the phase is 
already bigger than 27r, i.e. the fields are crossing 
the branch when going from one vacuum to the 
other. We can give an estimate of why this is 
the case. Given the boundary conditions we are 



setting for our field profiles, it is clear that, at the 
centre of the wall (i.e. z = 0) 

Po _ k Oio _ k _ 

27r ~ 2Nc' 27r ~ 2iVc 2 ' ^ ^ 

where, as usual, k = X^i"^*- Note that, by 
definition, the phase of the logarithm (given by 
{Nc — Nf)(3 + X^i'^i' modulo 27r) cancels here. 
It immediately follows that the evolution of the 
phases of the fields from z = towards the vac- 
uum, when TO is large, follows a very precise pat- 
ternQ: first, the phases of the matter condensates, 
ai change, more or less at the edge of the centre 
of the domain wall (defined by z — ±l/m); af- 
ter that, the gaugino condensate phase, /?, starts 
to change as well, to compensate the previous 
changes in such a way that the phase of the log- 
arithm cancels again at the vacuum. This can be 
formulated by saying that, when z = l/m, the 
different phases are given by 

— — ^ f^li - ^ _ ^ (21) 
2tt ~ 2Nc ' 2tt " 2Nc ' ^ ' 

■^This point is actually straightforward to realizeonce the 
large mass regime has been explained, see Eq. (p4). 
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This is the point at which the phase of the loga- 
rithm wiU acquire its biggest absolute value, be- 
fore it starts decreasing again. This is given by 
the changes produced both at ±l/m, i.e. 



I A (phase log)| ^ ^ , 



(22) 



Therefore, in order to have a change of the phase 
less or equal to 27r, = 1 (remember that k must 
be an integer). 

In practice it means that, for equal windings, 
only models with Nf = 1 will give domain walls 
profiles with no crossing of the branch. This im- 
plies a series of conceptual problems which are 
beyond the scope of this talk, but which should 
certainly be addressed at some point. 

Now that the main result has been presented, 
let us try to understand why the different situ- 
ations arise. In order to do that, we consider 
the cases where this model could be described by 
just one of the two condensates, i.e. either by the 
gaugino or by the matter condensate. This cor- 
responds to the limits of very small and very big 
mass, respectively. 

• For small masses (i.e. m <ti A) one can inte- 
grate out the gaugino condensate S and describe 
the theory in terms of the matter condensate QQ 
only. The resulting model accounts for the results 
we obtained for Nf = 1 jl^ and also for the up- 
per branch of the Nf — Nc — 1 case, studied by 
Smilga et al (lO). 

One can easily check that the second, lower 
branch present in this latter case corresponds to 
S ^ 0. If we have a look at the constraint equa- 
tion (Q9h evaluated at i?o the result is 



Po 



Nc f Nf 

cos TT 

Nf \ N, 



(23) 



The immediate conclusion one can draw is that 
there will only exist domain wall solutions with 
Po > if and only if ^ > 1/2. It has been argued 
in Ref. that the existence of this lower branch 
is directly connected to the fact that the Kahler 
metric for the S field is, in this case, singular. 
In that case, this lower branch would be nothing 
but an artifact due to the choice of metric. It 
should also be mentioned here that, according to 
Rcf. ||ll|], the total number of BPS-saturated do- 
main walls for a given value of Nf/Nc can only 



change by a multiple of two, when m is continu- 
ously varied, and that is the reason why the upper 
and lower branches end up collapsing into a non- 
BPS domain wall above m,. However this does 
not explain why they have to collapse at all. 

• For large masses (i.e. m ^ A) one can in- 
tegrate out the matter condensate (QQ) and our 
model will be described by the gaugino conden- 
sate only (i.e. we recover SUSY gluodynamics) . 
In order to do that, it is worth noticing that, in 
this large mass regime, M ^ S, and we can do 
the following identifications 



p{z)e 



ia.{z) 



p{z)e 



ioi.{z) 



i?(z)e*'^(^) , z<-l/m 



i?(z)e*(''(^)"2^) , z > 1/m 



(24) 



That is, to the left of the centre of the domain 
wall, both condensates behave in the same way 
whereas, to the right of it, their phase difference 
is 2ti. This is illustrated in Fig. 3, where we plot 
the phase difference between the two condensates 
(in units of 27r) as a function of the spatial coor- 
dinate, for several values of the mass parameter 
TO. One can also analyze the constraint equation 
( [l9| ) in this large mass limit, which will be given 
in terms of the modulus of the condensate i?o as 



i?o(l-ln(i?o)) 



Nc 



(25) 



A quick glance at this equation tells us that, in 
order to have i?o < 1 , which corresponds to finite- 
energy, well-defined domain walls, we must be in 
the case where ^ < 1/2. In other words, it does 
not seem possible to reach the to ^ oo limit, 
continuously from small to, in a model where 
^ > 1/2. The fact that, in that case, there 
are no BPS-saturated domain walls at large to, 
whereas we had two branches of solutions at small 
M explains why, at some point, those must have 
annihilated each other at to* . 

On the other hand, Eq. ( [25| ) is telling us that 
there is an analytic limit to — > cxd limit that we 
can construct when Nf/Nc < 1/2; using Eq. (p^ ) 
we can write down two BPS equations for S, one 
for the right hand side of the domain wall and 
another one for the left hand side. The solution 
is given by the thick line in Fig. 1 and, as we 
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Figure 3. 

Plot of the phase difference between condensates, (/3 — a) 
(in units of 2tt), as a function of the spatial coordinate, 
2, for Nc = 3, Nf = 1 and different values of the mass of 
the matter condensate: m = 2 (dotted), m = 20 (dash- 
dotted), m = 100 (dashed), m = 200 (soUd). 



Figure 4. 

Argand diagram of the profiles of the different condensates 
for Nc = 3, Nf = 2. Here wi = 1, L02 = 0, m\ = m2 = 
250. The gaugino condensate, 5, is given by the solid 
line while M\ and M2 correspond to the dot-dashed and 
dashed lines respectively. 



can see, coincides very well with what we would 
expect for the large m limit by looking at the 
finite mass results. 

4.2. Different windings 

Now that we have presented and explained 
most of the results in the literature, let us step 
onto new ground. As mentioned above, up to 
now the standard procedure was to consider all 
the matter fields transforming in the same way 
when going from one vacuum of the theory to 
another. We shall now break this degeneracy and 
assign different winding numbers to the different 
matter fields. At this stage we present results 
concerning a particular example, that of Nc = 3 
and Nf — 2, while a more detailed and general 
analysis will be presented elsewhere fl^ . 

According to the results of the previous subsec- 
tion, when A^c = 3 and Nf = 2, then A^//A^c > 
1/2 and we should expect two branches of BPS- 
saturated domain walls that annihilate each other 
at TO* (see Rcf. [0). And that is indeed what 
happens when uji = L02 = \. However, if we 



consider a different choice for the windings, i.e. 
cji = 1, clI2 = 0, then it turns out that one 
can construct BPS-saturated domain walls for 
any value of the mass to. The profiles of these 
fields can be seen in Fig. 4, where we plot the 
Argand diagram for S', Mi and M2 with masses 
TOi = TO2 = 250. There it can be seen that the 
matter condensate with ui — 1 follows a different 
path than the one with no winding, which sim- 
ply mimics the gaugino condensate, S*; and it is 
exactly the same path described by the matter 
field for Nc — Nf = 1. Actually it is possible 
to increase the value of TO2 while keeping toi fi- 
nite and constant, in order to integrate out the 
second condensate and recover the results we ob- 
tained for Nc — 3, Nf = 1. The key point to per- 
form this exercise is to notice that the relevant 
ratio that classifies the solutions is not Nf/Nc 
but k/Nc where, as defined before, k = J2i ^i- In 
the previous case, k = Nf, whereas now we have 
Nf = 2 with fc = 1. 

Therefore it should be in principle possible to 
construct BPS-saturated domain walls for any 
value of Nf, Nc, just by choosing the windings 
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of the matter fields in such a way that their sum, 
k, over the number of colours, Nc, is less than 
1/2. However it must be pointed out once again 
that, only when |fc| < 1, the paths described by 
the different fields between two vacua do not cross 
the logarithmic branch. Those are the only phys- 
ically meaningful domain walls that we know, so 
far, how to construct. 

5. Conclusions 

Let us summarize first the results obtained, 
before discussing what would be interesting to 
pursue in order to improve our understanding of 
these objects. We have looked in detail at domain 
walls in SUSY-QCD, as a way of understanding 
SUSY gluodynamics, the nearest relative to ordi- 
nary QCD. Using the TVY effective Lagrangian, 
we have constructed BPS-saturated domain walls 
for certain values of Nf, the number of matter 
condensates, and m, the mass of those conden- 
sates. The results so far obtained can be summa- 
rized as follows: 

• All matter fields transforming in the same 
way (equal windings) 

- If Nf/Nc < 1/2 there are BPS- 
saturated domain walls for any value 
of m. 

- If Nf/Nc > 1/2 there are BPS- 
saturated domain walls up to m». 

- The logarithmic branch of the scalar 
potential is not crossed only when 



N. 



f 



1. 



• Non-degenerate flavours (different wind- 
ings): any choice of coi such that k — 
^ ■ oJi = 1 gives BPS-saturated solutions for 
any {Nf,Nc) and m with no crossing of the 
branch. 

• The analysis of the constraint Im(e''''W) = 
const at the origin {z — 0) is the key to 
understand the results. 

So, at least, we have now a criterion to con- 
struct well-behaved BPS-saturated domain walls 
but, of course, this project is far from over. The 



next issue that should be addressed is that of the 
branches associated with the logarithm in the po- 
tential. As mentioned at the end of section 2, sev- 
eral suggestions have been put forward but none 
has so far provided us with a satisfactory answer. 

Another issue which is rather controversial is 
the dependence of these results on the choice 
of Kahler metric. In fact it has been widely 
claimed in the literature that the origin of the 
lower branch found in Refs ||l^ is a Kahler met- 
ric which blows up near the centre of those solu- 
tions. This is not a problem that one encounters 
in theories with higher supersymmetries, where 
the Kahler function is totally determined by the 
symmetries of the theory. Therefore an obvious 
way of trying to understand this dependence is 
through the connection of the N = 1 SUSY con- 
tructions here presented with those attempted 
with higher supersymmetries (for example, those 
of Refs. 0,|l|l). 

Finally it would be desirable to connect the re- 
sults obtained here with constructions of BPS- 
saturated domain walls in the large Nc limit in 
the context of SUSY gluodynamics, which have 
been performed in Ref. [Q. As mentioned in the 
introduction, these objects could play a very im- 
portant role in the D-brane description of SUSY- 
QCD. 
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